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Abstract 

Based on spherically symmetric reduction of loop quantum gravity, 
quantization of the portion interior to the horizon of a Reissner-Nordstrom 
black hole is studied. Classical phase space variables of all regions of 
such a black hole are calculated for the physical case > Q^. This 
calculation suggests a candidate for a classically unbounded function of 
which all divergent components of the curvature scalar are composed. The 
corresponding quantum operator is constructed and is shown explicitly 
to possess a bounded operator. Comparison of the obtained result with 
the one for the Swcharzschild case shows that the upper bound of the 
curvature operator of a charged black hole reduces to that of Schwarzschild 
at the limit Q — >■ 0. This local avoidance of singularity together with non- 
singular evolution equation indicates the role quantum geometry can play 
in treating classical singularity of such black holes. 

1 Introduction 

While describes large scale structures of the universe remarkably precise, Ein- 
stein's theory of general relativity (GR) fails to give a consistent picture at 
certain space-time points, called singularity, where physical observables of the 
theory, such as energy density, become divergent. It is widely believed that the 
true cure of such physically meaningless states of gravity predicted by GR, lies 
out of the domain of validity of the theory; it must be resolved in a regime in 
which both gravitational and quantum mechanical effects are dominant. 

Various attempts to construct a consistent theory describing such a regime 
has been made among them loop quantum gravity (LQG) [U [2] presents a clear 
picture of quantum space-time at Planck scales in the kinematical level. The 
program rests on canonical quantization of GR expressed as an SU(2) gauge 
theory on the 3 dimensional manifold S of arbitrary topology, where a topology 
of space-time M of the form M = Mx E is assumed, in a background independent 
manner. In such formulation of GR, the gravitational field is described by a pair 
of conjugate variables {A,E), where AJ,(a;|^is an SU{2) connection and E°;{x) 
is the densitised triad vector field, conjugate to A: 

-"^We are following Rovelli[l]'s notation in which indices fj,,u,... denote space-time, 
3-space and a,b, ... R'^indiccs. 



{A\{x), E'^ix')} = SnjGS^S'Jix - x'), (1.0.1) 

with G the gravitational constant and 7 the Immirzi parameter. 
The conjugate pair are constraint to satisfy the system 

= DaE!^ = 0, 
Hb = E°^F^j^ = 0, 
H - e,,kF:,E]El - 2(1 + -f')KmE!lE'^ = 0, 

which are called Gauss, spatial diffeomorphism and Hamiltonian constraints 
respectively. In fact, the task of finding a metric satisfying the Einstein's equa- 
tions, describing configuration of a gravitating system, is now replaced by finding 



a conjugate pair {A,E) satisfying the constraint system 1.0.2 

On quantization, one smears the basic fields (A, E) to holonomies of Af 

along a curve 7, defined by h^[A\ = P (^exp J^A^, and fluxes of E°' through 

the surface 5", defined by £^^(5") = Jg (PaUaE'^ . They form the holonomy-flux 
algebra in which holonomies act by multiplication, and fluxes act by derivation. 

Using a functional representation of quantum field theory and representing 
states as functionals of the cylindrical functions of holonomies, the kinematical 
Hilbert space of the theory is constructed. After imposing Gauss and diffeo- 
morphism constraints as operators on such states, the true gauge and diffeo- 
morphism invariant states of the theory turns out to be spin networks acted 
upon by holonomies and fluxes operators which form a unique representation. 
More precisely, a spin network is a triplet (F, j;, i„) consisting of a graph F with 
nodes in S, labeled by intertwiners «„, and links connecting different nodes, la- 
beled by SU{2) representations j/. The quantum geometrical picture suggested 
by LQG is manifest in quantization of geometrical observables, such as area 
and volume, as quantum operators acting on spin network states which result 
in discrete spectra and reflect the discrete nature of space-time. 

While the kinematical construction of the theory is well understood, the 
dynamics is still poorly under control. This can be investigated either within 
the spin foam formalism, a path integral analogue for time evolution of spin 
networks, or by studying solutions to the Hamiltonian constraint, which requires 
criteria for choosing between different variants and orderings of the operator, 
and consequently building the dynamical Hilbert space. 

Despite difficulties in solving the dynamics of the full theory, different at- 
tempts have been made in studying simpler situations which would shed light 
on the path toward complete understanding of the theory. To study specific 
gravitational configurations characterized by specific macroscopic spatial sym- 
metries, one should impose such symmetries suitably on the dynamical level of 
the quantum theory and find the true symmetric sector of the full theory. As for 
the Hydrogen atom in quantum mechanics, for instance, one would investigate 
the full quantum mechanical description of the system only by finding states 
satisfying the suitable dynamical equation, spherically symmetric Schrodinger 
equation, which form basis of the dynamical Hilbert space. However, in LQG 
due to difficulties in solving Hamiltonian constraint, one employs the symmetry 
reduction methods [Hj. In such models the desired symmetries are applied at 
classical level by choosing a symmetric pair { A, E) satisfying constraints or it 
can be imposed at the kinematical level of the quantum theory. Quantization 



of the symmetric pair will then be done which, evidently, gives rise to a simpler 
procedure than in full theory. This is in fact the idea of mini superspace which 
has been used extensively in quantum cosmology [4]. 

For the homogeneous and isotropic cosmological models, loop quantum cos- 
mology [5] has led to remarkable results such as replacing big bang singularity 
with a big bounce (Hj , modification of Fredmann equation [7 and incorporation 
of an inflation mechanism [5]. However, because of the reduction of degrees of 
freedom to a finite number, in fact one, the field theoretical aspects of the gravi- 
tation are lost; since all points of space are equal, spin network just has one node. 
The same method has been applied to the Schwarzschild black hole, since the 
interior region of such a black hole is a homogeneous space of Kantowski-Sachs 
type, and its singularity is shown to be resolved [9 . 

Intuitively, singularity resolution occurs as a result of fundamental discrete- 
ness of space; while in a classical continuum, divergences emerge as distance 
goes to zero, there is no room for divergences in quantum level since there is 
no zero distance below the Planck length. Put it slightly different, in the same 
manner quantum mechanical model for the hydrogen atom prevents a classical 
electron from collapsing to the nucleus, and hence energy divergence, by putting 
a lower bound on the spectrum of energy, quantum geometrical model put an 
upper bound on the classically divergent curvature components and prevents 
black holes from forming a singularity. 

In the spherically symmetric case, degrees of freedom of the theory cannot 
be reduced to a finite number. Instead, we are dealing with a midi-superspace 
having C/(l) as the gauge group. This makes calculations simpler than the full 
theory while still respects the field theoretical aspects of spherically symmetric 
quantum gravity. 

In this paper, we intend to study the quantum behavior of a Reissner- 
Nordstrom black hole. To this end, we employ spherically symmetric quantum 
geometry methods [101 E] to quantize the interior portion of such a space-time 
within which the classical singularity occurs. In section 2, we introduce the 
states and operators used in spherically symmetric quantum geometry. Section 
3 reviews the classical charged black hole metric, its horizons and singularity. 
And finally resolution of classical singularity will be discussed in section 4. 



2 Spherically Symmetric Quantum Gravity 

In this section, we outline the crucial points regarding the foundations and 
results of spherically symmetric quantum geometry developed by Bojowald [101 
[TT] . The reader is referred to the original articles for many technical subtleties 
which are omitted inevitably. 

Viewing LQG as a method to quantize connections, one would be able to 
impose a symmetry through two avenues: (i) to pick, in the classical level, 
only those connections which are invariant under symmetry group action and 
consequently reduce the phase space, and (ii) to restrict the distributional states 
of the quantum theory, at the kinematical level, only to invariant connection. 
It is shown in [TD] that both methods just mentioned lead to the same result 
in the case of spherical symmetry. Therefore, we proceed by describing the 
former-that of quantizing the classically reduced phase space. 



2.1 Classical Phase Space 



The manifold E having the symmetry group S with an isotropy subgroup F, 
can be decomposed asT, = T,/SxS/F. The connection can generally be written 
as A = A-^/s + As/p . Then A-^/g can be considered as the connection of the 
reduced theory and its liolonomics along curves in S/S' can be quantizcid. 

For the spherically symmetric case, S = M x 5^, and the symmetry group 
is 5 = SU{2). This implies identifying T,/S with M and the gauge group of the 
reduced theory F with U{1). Therefore, reduced connections are U{1) gauge 
fields on M. Roughly speaking, spherical symmetry implies that our basic fields, 
in the spherical coordinate (a;, 6, (j)), are independent of angular variables. Thus, 
the connection A(x) is just a function of the radial coordinate; A = A{x). 

The connections and triads of the reduced spherically symmetric phase space 
have the general form: 

A = Ax{x)T-idx + {Ai{x)Ti + A2{x)T2)d0 + {{Ai{x)T2-A2{x)Ti)sine+T3COSe)d(j), 

(2.1.1) 

E = E'={x)T3sin6d^ + {E^{x)ti + E'^{x)T2)sinede + {E^{x)t2 - E'^{x)Ti)d4„ 

(2.1.2) 

where Tj = — are the generators of su{2) algebra. They define the symplectic 
structure: 

O = ^ / dx{dA^ A dE"" + 2dAi A dE^ + 2dA2 AdE^). (2.1.3) 
27G J 

However, a suitable canonical transformation can be made resulting in canon- 
ical variables {A^{x),E'={x)), {-fK^{x),E't'{x)) and {r]{x), P'^ix)): 



n = [ dx{dA^ A dE"" + ^(7^0) A dE'l' + 2dr} A dP''), (2.1.4) 
27G J 

with being the (j) component of the extrinsic curvature. 

The Gauss constraint, generating [/(I) gauge transformations, takes the 
form: 

G[X] = j dxXiE"" + P") = 0, (2.1.5) 

where prime denotes differentiation with respect to x. 

Note that in terms of these variables, conjugate pair is not simply the 
connection-flux pair which suggests a different situation than the full theory. 

The Hamiltonian constraint can be written as 

^ ~2^ / '^''^^'^^^i^ (^^ ^ + ^^^^^^ + ^K^E^iA. + V') + 2i?^r; 

(2.1.6) 



2.2 Quantization 

Along the standard lines of constructing basic operators and states in the kine- 
matical Hilbert space of LQG, one starts with holonomies of the connections. 

Holonomies of A^ along curves 7 in K are defined as /i*^'*^ = exp ( |/ Ax{x)] 



which are elements in U{1). For point holonomies exp{inA^(x)) are used 
which belongs to the space of continuous almost periodic functions on the Bohr 
compactification of real line, and point holonomies of 77 € S^, have the form 
exp [irj{x)) which are elements of U{1). 

The kincmatical Hilbert space of the present reduced theory is the space 
spanned by spin network state Tg^k,^i 

(2.2.1) 

For a given graph g, these are cylindrical functions of holonomies along edges 
e of g. Vertices V{g) of such spin networks are labeled by irreducible K_Bo/ir 
representations //^ € M and irreducible representation fc,^ € Z, while edges 
are labeled by irreducible representations of U{1). 

Holonomies act on spin network states by multiplication. Their correspond- 
ing momenta, on the other hand, act by differentiation 

E^x)Tg^k..t. = 7y (fce+(.) + fc,-(,))r3,fe,^, (2.2.2) 
/ dxE'^{x)Tg^k,t. = ill ti.Tg^k^^, (2.2.3) 
I dxP^x)Tg^k,tj. = 27^2 ^ fc.T^^fc,^. (2.2.4) 



The volume operator can be expressed as = 47r J dx\E'>' {x)\^ \E'^ {x)\ 
which is diagonal in spin network representation 



Implementing the Gauss constraint as an operator on spin networks to select 
the gauge invariant states, leads to a restriction on labels 



G[A]T<,,fe,^ = 0=^K = - J(fce+(,) - fce-(.)). (2.2.7) 



Imposing this on 2.2.1 results in the gauge invariant states 



^9:fe,M = n e^P Q^ey"rfa;(Aa;(a::) + r;')^ Y[ exp{i^iulK^{v)) . (2.2.. 



£69 ^ ' y€V(g) 



3 Classical Reissner-Nordstrom Black Hole 



The Reissner-Nordstrom metric is a spherically symmetric solution to Einstein- 
Maxwell equations describing the space-time of a source with mass M and charge 
Q reads (in coordinate system (x, 9, 4>)): 

ds^ = - 1 + ^]df+{l + ^ dx^+ x^dn^. 3.0.9 

\ X / \ X / 

Horizons appear where goo — 0: 

x'^ - 2Mx + = — >x±=M± v/A/2 - Q2 (3.0.10) 

For the physically interesting case of > x^ , event horizons partition space- 
time into 3 regions: I (a; > x+), II (a;_ < x < x_|_), and III (0 < x < a;_). By 
inspecting the sign of goo, one observes that in region II, x and t interchange 
their roles and becomes time-like and space-like respectively. This defines a 
Kantowski-Sachs type homogeneous space. Region I and III, on the other hand, 
carry spherical symmetry. 



3.1 Classical Phase Space Variables 
3.1.1 Region II 

In region II, the metric of space-time takes the form 



= — -^-i) 'dt^+C^-^l.i^dx^ + fdn^ (3.1.1) 



According to definition of tetrad (frame) fields, 

5^1' = Vije^ei, (3.1.2) 
they can be determined only up to a Lorentz transformation. This leaves us 



with an 5*0(3,1) freedom in choosing tetrad. In fact, given the metric 3.1.2 



we are free to choose their sign and Minkowski indices, which can be viewed as 
sort of a labeling 4 tetrad fields. However, in order to serve as the fundamental 
fields for constructing the conjugate pair {A,E), a particular labeling must be 
chosen which will be clear below. 

The suitable choice for labeling 4 orthogonal frame fields turns out to be: 

e° = ±(— --^-Ij dt]e^ = ±tsined(j} 

= ±tde ; = ± - ^ - 1^^' dx, (3.1.3) 
which leads to compatible spin connection components: 

,30 _ , m _ f ™^ . , ,20 _ , ,02 /^2m ^ ^ 



a;- = -a;- = ( ^ - ^ Ux ; = = T " ^ " ^ ) '^^ (^.1.4) 



12 \ 1/2 

to — —LU 



,.01^1^^ ^ - 1 ) stnedcj) ; ^ _^2i ^ ^^^^^^^ ^3^^ 5) 



t 



The A field can be constructed using spin connections 

/2m -2 xi/2 



A'=Tl{f-'^) d., (3.1.6) 



A'=±j\^^-^-lj d0, (3.1.7) 
A' = ±j(—-^ - l) ' .sined(b, (3.1.8) 



A^ = ±cos0d(l3. (3.1.9) 

To construct the i? field on we choose a gauge in which e° = n^, the 
normal vector field to the spatial slice. This way we are in fact breaking the 
S0{3, 1) symmetry into 50(3) on a hypersurface with topology S = M x 5"^. 
The 3 triad fields become: 

/2to \1/2 
= ±tsined(p ; = ±tde ; = ± f — - - 1 j da;, (3.1.10) 



with determinant 



and inverse triad 



2 J2m Q2 \i/2 
det e = t^sirifii I — - - 1 



(3.1.11) 



1 1 /2m \-i/2 

= ^ - ± t ^« = = ± - ^ - V ^'-'-''^ 

The fields become: 

/2m 02 \^/2 
^2 = ±t ( — - ^ - 1 j 

= ±t^sined:,,. (3.1.13) 



The 3 triad fields 3.1.10 define their compatible spin connection, F'^ A + 
de^ = 

pi2 ^ _p2i ^ cos6'd0 (3.1.14) 

p3 ^ 1 (£312pl2 ^ ^321p21) ^ ^^..^^^^ (3_^_-^5) 

Extrinsic curvature is related to A via = ^4 — F 

Kl ^-Al^Ti^-'^^dx (3.1.16) 



7 " ' 



.2 _ 1 .2 , /^2m Q 



2 \ 1/2 



/^,^ = -A^ = ±(^— -^-1) (3.1.17) 



= ^ a; = ± - J - l) -ne (3.1.18) 

Note that had we chosen other Minkowski indices for tetrad [3". 1.31 we would not 

have obtained the conjugate pair(A, E) with correct indices satisfying {A]^{x), Ej{x')} — 

The phase space variables are determined up to a sign freedom. By demand- 
ing E and A to satisfy the diffeomorphism, Gauss and Hamiltonian constraints, 
their signs can be fixed relative to each other. All components of diffeomorphism 
and Gauss constraints are zero except: 

= 7^ - J - l) cosO {sgn{A\) + sgn{AlAl)] , (3.1.19) 
G2=d— -^-ij cose {sgn{El) + sgn{AlE^')] , (3.1.20) 



C = t' "" 



and Hamiltonian constraint gives: 

" ? " "'"'^ + sgn{Et)\ . (3.1.21) 

For the above constraints to be zero we must have: 

sgn(El) = -sgn{Et), (3.1.22) 

sgn{Al) = +1, (3.1.23) 

sgn{A\) = -sgn^Al). (3.1.24) 

This leaves us with two alternatives corresponding to the residual gauge freedom 
(6,P6) (-6, -ph). 

A A\ = cT^dr + bT2d9 + {cosOt^ - bsin9Ti)d(j) ^ ^c) 

\ Ef = pcT^sinddr + PbT2sindde ~ PbTid^, 

f Al^ ^ cT^dr - bT2d9 + {cos9t3 + bsin9Ti)d4> (3 126) 
\ E!^ ^ pcTssinddr - PbT2sin6de + pbTid^, 

where, 

Pc = ±i^P.-t(^- J-l)'^'. (3.1.27) 

The momentum pc ~ ztt'^ is a monotonic function and can be interpreted as an 
internal time parameter (as is interpreted in |5] for the case of the Kantowski- 
Sachs minisuperspace of Schwarzshild black hole). 



3.1.2 Region I and III 



The analogues calculations for regions I and III with line element 3.0.9 leads to 
the folowing phase space coordinates 



I 



II 



Al^ = cT^dr + hT2dO + {cosOt^ — bsin9Ti)d(j) 
Ef ^ pcT3sin9dr + PbT2sin9de - Pbnd^, 

A]^ = cT^dr — bT2d9 + {cos9t3 + bsin9Ti)d(j) 
Ef = pcT3sin9dr - PbT2sin9de + PbTid^, 



where, 



±71 



2ra 

X 



1/2 



T7 



Pc 



±x ; pb^ X [1 



2m Q 

X X 

This defines variables introduces in section \2A\ as 



TO 

X 

2\ 1/2 



A^ ^ d , E'-" ^ Pc ; 

jK^ = b , E'^ = Pb ; 
ri = (2n + l)7r , P** = 0, 



(3.1.28) 
(3.1.29) 



(3.1.30) 

(3.1.31) 

(3.1.32) 
(3.1.33) 



which constitute a 4 dimensional phase space. 

The singularity occurs at x = 0, or equivalently E^ — 0, where the scalar 
curvature become infinite and lies in the spherically symmetric region III. 



4 Local Singularity Resolution: Curvature Bound- 
edness 

In classical general relativity, a singularity is characterized by two criteria: (i) 
locally: the scalar curvature, made out of Riemann tensor, at a point of space- 
time becomes divergent, and (ii) globally: Einstein's equations cannot be solved 
for all values of t. Below, we consider the local criteria in the quantum theory 
and make a comment on the global one in the conclusion section. 

A look at Kretschmann scalar curvature of the Reissner-Nordstrom black 
hole [12 

reveals that the irremovable curvature singularity of a charged black hole oc- 
curs at X = 0. It lies in the spherically symmetric region III of the Reissner- 
Nordstrom space-time. Our goal is to show that the quantum operator corre- 
sponding to ^, on the other hand, exhibits a spectrum which is bounded above. 
To construct such an operator, we need the classical function to be expressed 
suitably in terms of well-defined operators in the reduced theory. '^s"-{e ) _ i 

is a candidate to serve as the desired classical function. However, one cannot 
naively replace the inverse squired of triad with its operator analogue since, as 
can be seen in 2.2.2[ it has zero eigenvalue. Nevertheless, being expressed as 



a Poisson bracket of functions having wcll-dcfincd operators, its corresponding 
quantum operator can be realized by replacing Poisson bracket with {ih times) 
commutator. This can be done by the simple relation: 



Since E'^{x) is diagonal on spin network states, ^J {x)\ can be defined as 
a diagonal operator whose eigenvalues are square roots of (absolute value of) 
those of E^{x). 

Along the inhomogeneous direction x, we can expand holonomy hx{x) as: 

hx = exp (^i J dxAx{x)^ « 1 + ieA^iv), (4.0.36) 

where e = — v is the coordinate distance between two vertices v and 1/+ 
connected by the edge e. This enables us to, as is usual in LQG, express TZ to 
the order of as: 



and consequently its corresponding operator as: 

1 



tr T3hx 



(4.0.37) 



(4.0.38) 



To calculate the above expression, we write holonomy as: 

hx = exp (^i j^dxAxix)^ = cos(^ j A^) + 2T3sin{^ j A^), (4.0.39) 
and therefore, 

K \h-\^J\E-{x)\ = ^\E-{x)\-cos{\ j Ax)^J\E-{x)\cos{\ j A,) 

-sm(^ j Ax)^J\E-{x)\sin{^ j A^) 
+2t3 {cos{^ j Ax)^\E^{x)\sin{^ J A^) 
-sm(i / Ax)^J\E-{x)\cos{^ ( A,) 



This leads to 
1 



n = 



2777^1 



PI 



cos(^ j Ax)^\E-{x)\sin{^ j A,) - sin{^ j A^)sJ\¥- 



(.t)|cos(- 
(4.0.40) 
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Figure 1: The upper bounded spectrum of ^ for a charged black hole 

To investigate its action on spin network states, one can use identities ex- 
pressing sinx and cosx in terms of exponentials and write: 



e.u 



which gives rise to the spectrum 



T^Tg^k^ti = 2-K^^Pl [v^l'^'=+(^) +ke-{x) +2| - \j]^\ke+{x) + K-{x) -2|j Tg,fc^^. 



(4.0.41) 



Such an operator exposes two features: (i) it is diagonal on spin network 
states and therefore commutes with all operators in kinematical Hilbert space, 
and (ii) it has a bounded spectrum (see fi gure 1) with maximum value of 
r- ^ „ — . Thus, the scalar curvature 4.0.34 which is classically divergent, at 

quantum level has a maximum value oii 



/ ASAPx^ - 96MQ'^x + 56Q2 

\ / max 

- -Air U'' - + ^VoA2) 

. This reflects the fact that the divergence of ^ is kinematically well-behaved 
at quantum level. 



5 Comparison with Schwarzschild Case 



The Kretschmann scalar curvature of the Schwarzschild black hole, W^^f"^ R^^p^ 

48M" 



t^, indicates that in this case the irremovable singularity of space-time lies 
at t = as well. However, the region within the event horizon of such a space- 
time is of a Kantowski-Sachs type and methods introduced in section [2] are not 
applicable to it. Loop quantization of such space-times, which are of interest 
in cosmology as well, is discussed in [9]. However, in [9l the issue of resolu- 
tion of classical singularity in quantum regime is addressed by illustrating the 
stronger condition of non-singular behavior of the quantum dynamical equation, 
and the boundedness of scalar curvature operator is not exhibited directly. Be- 
low, we sketch the construction of such an operator in a manner similar to the 
calculations of the previous sections. 

The classical phase space variables would be the same calculated in subsec- 
tion 



3.1.1 in the limit (5 = 0. Therefore the set {c,Pc} used in this section is 
meant to be understood as the Q = limit of [3.1.27] and |3.1.27| 
Consider the following quantity on the classical phase space 

Following the methods presented in [5] , expand the holonomy along x direc- 
tion of S = M X 5^ with oriented length r as 

= 1 e I dxcT3 + 0(6^), (5.0.44) 
Jo 



and write TZ as 



n = ^tr (ra/^M {h^J^-\ /h}) . (5.0.45) 



Now, quantization would be straightforward 



7^ = 



('^°*(y)^^i^^*"'^y) ^ s«'^(y)v^cos(y)) . 



(5.0.46) 



Its action on |t, /i) , which arc the simplified version of the spin network states 
in this reduced model (with /i being the oriented length along the equator of 
S"^), then becomes: 

n\T,^i) = —^{./W+^-^JW^\)\T,^i). (5.0.47) 

Such an operator has a bounded spectrum with maximum value of r- — . 

This means that the scalar curvature, which is classically divergent, at quantum 
level has a maximum value of: 



/48M2 



UP 



(5.0.48) 



Note that in the Umit Q — > 0, the maximum value of the curvature of a 
Reissner-Nordstrom black hole, 4.0.42[ reduces to that of Schwarzschild black 
hole 15.0.481 



6 Conclusion 



The remarkable result in quantum geometry is the discreteness of space mani- 
fested in discrete spectrum of geometrical observables such as area and volume 
operators. This enables us to express other classically ill-defined functions in 
terms of basic quantized fields, and to study their quantum mechanical behav- 
ior. One of such functions is the scalar curvature of a solution of Einstein's 
equations. In cosmological models, divergence of the scalar curvature is due to 
unboundedness of inverse scale factor at initial singularity, while in black holes 
this happens at the center of black hole where the function - , constituting all 
curvature components, becomes arbitrary large. In the above sections, we have 
constructed the quantum analogue of ^ for the spherically symmetric region 
interior to the horizons of a Reissner-Nordstrom black hole, and have shown 
that it is an operator possessing a bounded spectrum. However, this is just the 
local criterion for absence of singularity; the global criterion requires the quan- 
tum Einstein's equation representing the dynamics to be non-singular as well. 
To study the dynamics, we must encode time evolution coded in Hamiltonian 
constraint. To this end, one must assign a time parameter to the theory. Such 
a quest is usually realized via coupling a scalar field to the theory. However, 
in the case of spherical symmetric space of a charged black hole, we can con- 
sider the monotonic function E^{x) = ix^ as an "internal time". Since for 
the interior portion of a Reissner-Nordstrom space (region II) the internal time 
is taken as ±t^, this strongly suggests, for the sake of continuity, to take the 
same function ±a;^(= ±t^) to serve as this purpose. In the quantum level, thus, 
the eigenvalues of E^{x) play the role of discrete internal time. To see whether 
or not the time evolution of quantum space-time come to a halt at singularity 
E^{x) = 0, we must act by Hamiltonian constraint on spin network states and 
inspect whether or not the coefficients of its eigenvalues remain non-vanishing 
at such a point. The main task of realizing the Hamiltonian operator is done 
and extensively discussed in TT , and as pointed out in [13! , among the various 
options to choose an ordering for Hamiltonian constraint, the symmetric one, 
H = Hi + H,y , turns out to be the one with non- vanishing coefficients and can 
serve as the criterion of non-singular evolution equation. Therefore, absence of 
singularity can also play the role of a criterion for selecting the right ordering 
for the theory. 
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